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Abstract 

A phenomenological method based on the natural orbital representation is applied to construct 
the ground state one-body density matrix which describes correctly both density and momentum 
distributions in '''He, and ^"Ca nuclei. The parameters of the matrix are fixed by a best fit to 
the experimental density distribution and to the correlated nucleon momentum distribution. The 
method allows the natural orbitals, the occupation probabilities and the depletion of the Fermi sea to 
be obtained. Ground-state characteristics of ''ffe, ^^O and ^"Ca nuclei, such as rms radii and mean 
kinetic energies are calculated, as well. 



1 Introduction 

The one-body density matrix (OBDM) plays an important role in the nuclear structure theory. It is known 
that the nuclear wave function contains much more information than it is accessible to observation ^ . 
All measurable characteristics of the nuclear ground state (apart from the energy) correspond to single- 
particle operators whose expectation values can be expressed in terms of the OBDM elements. The 
ground state energy can be calculated by a sum rule analysis |[ [s] also by means of the OBDM. Thus 
the determination of the OBDM from first principles ||^ is one of the most essential aims of the nuclear 
theory. 
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As shown in |0, ^ (and discussed in |l^, |T^, |T^), an appropriate criterion for the proximity of the 
OBDM po(r, r') corresponding to a single-Slater determinant wave function (i.e. Pq=pq, for instance in 
the Hartree-Fock approximation) to the OBDM p(r, r') of the true (correlated) ground state is that the 
"mean-square deviation per particle" 

a ^ A-'Tr[{p - pof] (1) 

(where A is the mass number) should be minimal. It has been shown in |^ that a is minimal when po{r, r') 
corresponds to single Slater determinant wave function constructed with natural orbitals, i.e. with the 
single-particle wave functions which diagonalize the OBDM p{r,r') As was shown in [||, if in the 

Hartree-Fock approximation the OBDM diagonal elements pQ^{r,r') are fitted to reproduce correctly 
the exact density distribution p(r) = p{r, r) (i.e. pQ^{r, r) 2± p(r)), the extremum property (T=amin leads 
inevitably to an increasing deviation between the non-diagonal elements of these two matrices (p(r,r') 
and pQ^{r,r') at r 7^ r'). Th non-diagonal elements are related, however, to the nucleon momentum 
distribution (NMD): 

n(k) = yp(r,r')exp[ik(r-r')]drdr'. (2) 

This leads to the conclusion that generally the mean-field approximation is unable to give simultane- 
ously a correct description of the two basic nuclear ground state characteristics, namely the density and 
momentum distributions. This has been supported by the calculations in the case of ^He nucleus fl^ . 

The role of the nucleon-nucleon correlations in the nuclear systems can be seen by the analyses of 
the quantity a [lO[ |. While in the Hartree-Fock approximation a vanishes and the contribution of the 
random-phase approximation to a is about 0.002, for realistic N-N forces the short- and medium-range 
correlation effects lead to a value oiamin — 0.02-^0.03 that implies a Fermi sea depletion of about 10-^15%. 
The latter is confirmed by the experimental data The short-range and tensor correlations are 

responsible for the existence of high-momentum components in the realistic momentum distributions in 
nuclei which is not the case in the mean-field approximation. 

A phenomenological method to construct a more realistic OBDM is suggested in [|l], |2) for the cases of 
^^O and ^'^Ca nuclei, both with fractional occupation numbers and full occupation of the first A natural 
orbitals. The latter are expanded in terms of harmonic-oscillator wave functions. The expansion coeffi- 
cients, as well as the occupation numbers are fixed by a best fit to the experimental density distributions 
of ^^O and '^'^Ca minimizing the rms percentage deviation between the experimental and trial radial 
density. A satisfactory description of g.s. properties, such as the g.s. energy and rms radii, is obtained 
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using the OBDM constructed within this method, in contrast with the predictions of the Hartree-Fock 
approximation. 

Here we would like to emphasize that according to the discussion given above, the OBDM obtained 
in the method from |0, |j cannot be a realistic one because it is obtained by a best fit only of its diagonal 
elements to the experimental density distribution (which is not very sensitive to the N-N correlations) 
and is not aimed to describe the nucleon momentum distribution which is related also to the non-diagonal 
elements of the OBDM and is mostly affected by the presence of short-range and tensor N-N correlations 
in the nuclear system. 

The aim of this work is to construct the g.s. OBDM in '^He, ^^O and ^"Ca nuclei in a consistent way 
following the method from ^ but providing in an optimal way a correct description of both density 
and momentum distributions in nuclei considered. Solving this problem, we show that, in principle, the 
expansion of the natural orbitals in terms of harmonic oscillator functions (in the truncated s.p. space) 
cannot give a correct description of the realistic high-momentum components of the nucleon momentum 
distribution. This imposes to use in the expansion of the natural orbitals another set of s.p. functions, 
for instance that one corresponding to infinite square-well potential. 

The method suggested in this work enables us to construct a more realistic OBDM which includes 
that part of the nucleon correlations which are responsible for the high-momentum components of the 
nucleon momentum distribution and is consistent with the nuclear density distribution. 

Some basic relations are given in Section 2. The phenomenological expression for the OBDM is 
presented in Section 3. The results of the numerical calculations are given and discussed in Section 4. 

2 Nuclear ground-state properties 

In order to evaluate the single-particle characteristics of the nuclear ground state one needs the OBDM 
which is defined by the expression: 



where 4'(ri, r2, r^i) is the normalized A- nucleon ground-state wave function. The integration in Eq.(3) 
is carried out over the radius vectors and summation over spin and isospin variables is implied. 

The one-body density matrix has the simplest form in the so called natural orbital representation p3| : 




(3) 




(4) 
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where Xi are the occupation numbers (eigenvalues of p) corresponding to the natural orbitals ipi (r) (eigen- 
functions of p) which form a complete orthonormal set. The values of satisfy the general conditions: 
< Ai < 1 and Xi — A. Usually there are A orbitals ipi (r) for which the occupation probabilities 

i 

Xi are significantly larger than those for the others. As in the mean field approach, these are called the 
hole-state orbitals while the others are called particle-state orbitals p7| . 
In the momentum space the OBDM (Eq.(3)) reads 

n{k,k')=A J ^'*(k,k2,...,kA)*(k',k2,...,k^)dk2...dk^, (5) 

where ^(ki, k2, k^) is the Fourier transform of the A-nucleon wave function ^'(ri, r2, r^i). Its 
diagonal elements determine the nucleon momentum distribution (Eq.(2)). The associated natural orbital 
representation 

n(k,k') =5]A,Vi*(k)V;.(k') (6) 

i 

depends on the same occupation numbers Xi entering Eq.(4) and the Fourier transform of the natural 
orbitals ipiir) in the momentum space 



(k) = (27r)-3/2 j Vz(r) exp(ikr)dr. (7) 



For the purposes of our work a particular attention will be paid to the local density distribution (the 
diagonal elements of the OBDM in the coordinate space): 

p(r)^p(r,r)=^A.|^,(r)p (8) 

i 

and the nucleon momentum distribution (the diagonal elements of the OBDM in the momentum space): 

n(k) =n(k,k) =^A,|^,;(k)|2. (9) 

i 

The second moments of p{r) and n(k) define respectively the rms radius of the nucleus 

{ry/' = f^A-' J pirydvY' (10) 



and its mean kinetic energy 



(T)^^ I n{k)edk. (11) 



The knowledge of a realistic one-body density matrix would allow to describe correctly the nuclear 
characteristics, such as the nucleon momentum and density distributions, rms radii and mean kinetic 
energies taking into account the effects of the nucleon-nucleon correlations in the nuclear system. 
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3 Phenomenological method for the constructing of the one- 
body density matrix 

In this Section a phenomenological expression for the OBDM is given. The parameters in this expression 
can be obtained by a numerical procedure by fitting the local density p(r) and the momentum distribution 
n{k) to the corresponding experimental data (or to realistic theoretical estimations) [|l^, |l^, 

For closed-shell nuclei with equal numbers of protons and neutrons it is reasonable to assume that the 
protons and neutrons have the same one-body density matrix. For such nuclei (with total spin J = 0) 
the OBDM has to be diagonalized in the {Ijm} subspace of the complete space of the natural orbitals 

ill (i 

,j,m being the quantum numbers corresponding to the s.p. orbital and total momentum and its 
projection). Due to the spherical symmetry of the considered nuclei the natural orbitals can be looked 
for in the form: 

(r) = V„;m (r) - Rni (r)!;™ {0,ip). (12) 

It is known that the radial part of the natural orbitals i?„/ (r) and especially of the particle-state orbitals 
differ from the single-particle wave functions obtained within the mean field approximation to the ground 
state of a fermion system |l^, |l^, ^ . Following exactly the method suggested in jl[ || , we expand 
the radial part of the natural orbitals of the OBDM in terms of three single-particle functions preserving 
all the usual symmetries for spherical nuclei: 

3 

Ra{r) = J2G:Mr), (a = nl) (13) 

1=1 

In Eq.(13) {fni{r)} is a set of orthonormal s.p. wave functions. The expansion coefficients Cf are fitting 
parameters satisfying the orthonormalization conditions. For the numerical calculations it is convenient 
to re-express them using polar coordinates: 

Ci" = cos 0a 

C2 = Sm0a cos ipa (14) 

C3 = sin 6*0 sin (/3a . 

In this way the normalization is automatically guaranteed. The orthogonality among states of the same 
symmetry reduces the number of parameters involved in Eq.(14). 

By means of Eqs. (12-14) one can obtain the expression for the one-body density matrix (Eq.(4)) 
within this method. The diagonal elements of the OBDM in coordinate space determine the local density 
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distribution: 

'^W = iE4(2^ + l)A„iKi(r-)p. (15) 

nl 

The occupation numbers A^^ satisfy the condition: 

^4(2; + l)A„, =^ (16) 

nl 

and the normahzation of p{r) is: 

I p{v)dv - A. (17) 
An expression similar to that for p(r) determines the nucleon momentum distribution: 

n(fc) = ^5l4(2/ + l)A„i|-R«i(fc)l'; Jn{k)dk^A. (18) 

nl 

In Eq.(18) 

3 

Ra{k) = Y,C7Mk), {a = nl) (19) 

1=1 

where (pni{k) are the orthonormaHzed s.p. wave functions in the momentum space. Obviously, due to 
the truncation of the expansion (13) (and the corresponding one (19)) up to the radial part of the 
natural orbitals depends strongly on the choice of the single-particle wave functions {ipni}- We use three 
sets of s.p. wave functions {^ni\ corresponding to: 1) the harmonic-oscillator potential (HO): 

V[r) = -Vo + ^mw^r^ {Vq > 0), (20) 

2) the square-well potential with infinite walls (SW): 

f -Vo, r<x, Vq>Q 
V{r) = <^ (21) 

I oo, r > X 

and 3) the modified harmonic-oscillator potential (MHO) pst : 

V{r) = -Vo + \mv?r^ + Vo>0,B>0. (22) 
2 

The modified harmonic-oscillator potential (Eq.(22)) behaves like a harmonic-oscillator one for large 
values of r but it has in addition a repulsive term {B determines the strength of the repulsion) which 
is the dominant one at short distances from the origin. One can expect that this term simulates to 
some extent the effects of strong repulsion in the nucleon-nucleon interaction at small distances. The 
single-particle potential (22) had been suggested and used |2^, ^ for the study of charge formfactors 
and nucleon momentum distributions of light nuclei. One of the main advantages of this potential is that 
analytical expressions are derived for the single-particle wave functions and for other useful quantities 



(see e.g., |^^). Though the MHO potential (22) leads to an infinite repulsion at the origin (r = 0), the 
corresponding single-particle wave functions which are used in the present work are correct for all values 
of r. 

In the natural orbital representation the local density and the momentum distributions (Eqs.(15) and 
(18)) are expressed as sums over the whole single-particle space. In the numerical calculations we include 
in these sums all hole-state orbitals and several particle-state orbitals. In this way the expression for the 
OBDM is determined and it can be used within a minimization procedure which leads to the constructing 
of a realistic one-body density matrix. 



4 Numerical calculations. Results and discussion 

The phenomenological method for constructing the realistic one-body density matrix given in the previous 
Section is applied to the case of the double-closed shell nuclei ^He, ^^O and ^°Ca. Obviously, the 
expansion (4) should be truncated. In the case of '^He and ^^O we include the Is, Ip, Id and 2s 
states and for '^^Ca nucleus-the Is, Ip, Id, 2s, 2p, 1/ and \g states. The parameters of the OBDM, 
namely the occupation numbers {A„;}, the coefficients {C,"'} and the parameters of the sets of the 
orthonormal single-particle functions corresponding to the potentials (20)- (22) (the oscillator parameter 
a = [mw/h)^/'^ , the radius of the well x and the parameter b — SmB/fi^) are fitting parameters. The 
values of the parameters are fixed by a best fit of the diagonal elements of the OBDM in the coordinate 
space and in the momentum space (Eqs.(15) and (18)) to the experimental local density p{r) obtained 



in |26, |2^, ^ and to the nucleon momentum distribution n{k). In the case of '^He we use the 
experimental data for n{k) obtained by means of the y-scaling analysis relying on the assumption 
that the 1/q expansion is valid. Theoretical estimations for n{k) in ^^O and '*"Ca obtained within the 
Jastrow correlation method (JCM) ||l^, |l^, ^ are used instead of experimental ones due to the lack of 
empirical data for these nuclei. Though the low-order approximation to the JCM has some deficiences, 
we use it due to the obtained simple analytical expressions for the one-body nuclear characteristics, in 
particular, for the n{k). The latter reproduce the results from the complicated exp(S)-calculations for ^^O 
[0 including the high-momentum tail of n{k). It is possible also more realistic theoretical predictions 
for n{k) in the nuclei considered to be used. 

Our fitting criterion is to minimize the relative deviation between the experimental and theoretical 
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local density distributions 



(23) 



where q is the number of points under which minimization is performed, and between the theoretical 
nucleon momentum distribution for ^^O and ^^Ca (or experimental n{k) for ^iJe) and that corresponding 
to the OBDM proposed in this work. As far as the latter is concerned it is preferable to consider Ign(fc) 
rather than n{k) because the values of n{k) rapidly decrease by several orders of magnitude with the 
increase of the momentum k: 



Otherwise, our minimization criterion would not be sensitive to both parts (low-momentum and high- 
momentum region) of the nucleon momentum distribution. Taking into account Eqs.(23) and (24) we 
minimize the expression: 



The minimization is accomplished using the program MINUIT. There have been involved in it 12 (or 
13 when MHO s.p. wave functions are used) fitting parameters in the case of '^He and ^^O and 20 
parameters for ^°Co nucleus. 

The quality of the minimization procedure is illustrated in Table 1. There are listed the values of Xd 
(Eq.(23)), Xm (Eq-(24)) per fitting point q and F (Eq.(25)). For all nuclei examined it is seen that the 
minimal value of F is achieved using square-well s.p. wave functions. For example, in the case of ^^O 
the local density distribution (15) and the nucleon momentum distribution (18) are approximated with 
almost the same accuracy. A low value of F is also obtained when modified harmonic-oscillator s.p. wave 
functions are used but Xd/q and Xm/q differ by two orders of magnitude. It means that in this case the 
local density distribution is described better then the nucleon momentum distribution. As is well known, 
the harmonic-oscillator s.p. wave functions are appropriate to describe the local density distribution of 
^^O within the simple shell model. Considering the values of Xd, Xm ^ ^^^^ proposed 

method in this case is unable to reproduce simultaneously the local density and the nucleon momentum 
distribution. Similar calculations can be drawn considering the minimization procedure for '^^Ca. The 
best fit value of F is obtained using square-well s.p. wave functions. The local density distribution and 
the nucleon momentum distribution are described satisfactorily to almost the same accuracy. 

The values of the occupation numbers Xni, of the coefficients Cf' and of the parameters of the s.p. 
wave functions in the expansion of the OBDM for ^He, ^^O and *°Ca are given in Tables 2, 3 and 4, 




-I 2 



(24) 




(25) 



8 



respectively. Thus the one-body density matrix is fully determined within the method proposed. The 
values of the occupation numbers A„; are compared with the occupation numbers r]ni p9[ obtained within 
the natural orbital representation of the OBDM in the low-order approximation of the Jastrow correlation 
method which can be considered as an "exact limit" of the method proposed in the present paper. It is 
seen that the values of A„; in the case of square-well s.p. wave functions are in good agreement with the 
natural occupation numbers r]ni. 

The depletion of the Fermi sea in *He is 18.6%, 13.2% and 11.3% in the cases when harmonic- 
oscillator, modified harmonic-oscillator- and square-well s.p. wave functions are used, respectively, in the 
expansions (13) and (19). Its corresponding values using the same types of s.p. wave functions in the case 
of are 18%, 16.6% and 12.4%, respectively. In the case of "^Ca the depletion is 20% and 15% when 
harmonic-oscillator and square-well s.p. wave functions are used, respectively. It can be seen again that 
the best description of the OBDM within the proposed procedure is achieved when square-well s.p. wave 
functions are used because the depletion in this case is minimal and it is the closest to the corresponding 
value in the natural orbital representation jl^]-6% for "^He, 4% for ^^O and 6.3% for ^^Ca and to the 
experimental one |^-9.4% for ''°Ca. We would like to mention here that the differences of our results 
from those obtained in [ p^ and from the experimental ones for ^^Ca ]3^ have to be considered bearing in 
mind the deficiences of the low-order approximation to the Jastrow correlation method used in , the 
model dependence of the analyses of the empirical data and, of course, the limitations of the suggested 
method. Nevertheless, we emphasize that our results are within the theoretical limits for the depletion 
of the Fermi sea which is expected to be approximately between 5 and 15%. 

In Figs.l and 2 the results for the local density and momentum distribution in ^_ffe, ^^O and '^'^Ca 
obtained by means of the one-body density matrix determined in the model are presented. The best 
solution is reached by using of square-well s.p. wave functions. As is seen from Fig.l an acceptable 
quantitative agreement with the experimental data for the density distribution p{r) (at least in the surface 
region) can be achieved almost independently on the choice of the s.p. wave function sets. However, only 
the use of the SW s.p. wave functions corresponding to the potential (21) gives a correct description of the 
realistic high-momentum components of the nucleon momentum distribution (see Fig. 2). As is expected, 
the use of HO s.p. wave functions does not reproduce realistically n{k). It can be also seen from Figs. 2b 
and 2c that although the OBDM obtained in the method from 0, || gives a correct description of the 
density distribution, it is unable to describe the high-momentum behaviour of the momentum distribution 
which is sensitive to the short-range N-N correlations. In the cases of ^He and ^^O (Figs. 2a and 2b) 
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the use of MHO s.p. wave functions corresponding to the potential (22) leads to an improvement of the 
behaviour of n{k) at higher momenta which is due to the existence of the repulsive term in the potential. 

In Table 5 the calculated rms r adii (r2)i/2 (Eq.(lO)) and mean kinetic energies per nucleon (T) /A 
(Eq.(ll)) of ^He, ^^O and "^^Ca and their values deduced from the JCM are presented. The values 
of the rms radii for all nuclei examined are in accord with the corresponding experimental values ||2^, ^ . 
Due to the fact mentioned above that the expansion of the natural orbitals in terms of HO functions 
cannot give a correct description of the realistic high-momentum components of the nucleon momentum 
distribution (see Fig. 2c), it can be seen in Table 5 a quite increased value of (T) /A for ^"Ca in comparison 
with the result obtained by using SW s.p. wave function set. Therefore, we have demonstrated that the 
SW s.p. wave functions are the best ones among the other used in the sense of correct description of 
both density and momentum distributions in '^He, ^^O and ^^Ca. Here we would like to mention that 
our choice of the truncation in the expansions (13) and (19) up to i = 3 was imposed from one side by 
the desire to follow exactly the method previously suggested in |^ in order to compare our results with 
those from the mentioned works. On the other hand, it is well-known that the expansions (13) and (19) 
in terms of the complete set of s.p. wave functions would be exact expressions for the natural orbitals. Of 
course, this is impossible to be done and in the practical applications we have to find a proper value of i 
which would lead to a reasonable number of the minimization parameters. The inclusion of higher terms 
in the expansions (i > 3) would make the results of the method less sensitive to the particular choice of 
the s.p. wave function set, but at the same time this would increase enormously the difficulties of the 
practical minimization procedure. Hence, the applications of the method would be hardly possible. In 
our opinion, the choice i = 3 is an acceptable compromise within the method. 



5 Conclusions 

The results of this work can be summarized as follows: 

i) a phenomenological method for a consistent construction of the one-body density matrix is suggested. 
In contrast to the method from ^ , our procedure allows to obtain an one-body density matrix which 
gives realistically both the local density and the nucleon momentum distribution of "^He, ^^O and ''"Co 
nuclei. 

ii) The method gives the natural orbitals, the occupation probabilities and the depletion of the Fermi sea 
in the nuclei considered. 

iii) The rms radii and mean kinetic energies of *_ffe, ^^O and ^°Ca are calculated with the method. 
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iv) Consideration of both non-diagonal and diagonal elements of p(r, r') makes it possible to treat in a 
consistent way the effects of the N-N short-range correlations in nuclei on the one-body density matrix 
and the related ground state characteristics. 
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Figure Captions 

Figure 1. Local density distribution in '^He (a), ^^O (b) and ^°Ca (c) obtained by using different sets of 
s.p. wave functions in Eq.(15). The normalization is J p{r)dr = A/2. 

Figure 2. Nucleon momentum distribution in ^i7e (a), (b) and ^°Ca (c) obtained by using different 
sets of s.p. wave functions in Eq.(18). The normalization is / n(k)dk = 1. 
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Table 1. x^-values per fitting point for the local density (Eq.(23)) and the nucleon momentum distribution 
(Eq.(24)) and the F-value (Eq.(25)) for ^He, ^^O and ^^Ca. 



Nuclei 


S.p. wave 




Xm/q 


F 




functions 










HO 


0.0040 


0.0379 


1.204 




MHO 


0.0031 


0.0260 


1.027 




SW 


0.0068 


0.0036 


0.932 




HO 


0.0098 


1.2458 


6.239 




MHO 


0.0019 


0.1271 


2.104 




SW 


0.0199 


0.0170 


1.888 




HO 


0.1805 


0.6678 


8.147 




SW 


0.0168 


0.0352 


2.190 
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Table 2. Occupation numbers Xni, expansion coefRcients C"' (i=l,2,3), parameters of the s.p. wave 
functions and natural occupation numbers rini p9[ for ^He. 

HO s.p. wave functions MHO s.p. wave functions SW s.p. wave functions 



nl a = 0.648 fm^^ a = 0.687 fm^i b = 0.511 x = 3.70 fm ijni 



0.9806 0.9699 0.9755 

Is 0.8138 0.1937 0.8679 0.2435 0.8873 -0.2191 0.9428 

-0.0308 0.0001 0.0182 

-0.5893 -0.9643 0.1411 

Ip 0.0392 -0.5011 0.0145 -0.1136 0.0063 -0.9576 0.0097 

-0.6337 0.2392 0.2513 

0.1848 -0.3959 -0.3405 

Id 0.0033 0.1256 0.0064 0.1115 0.0187 0.5185 0.0027 

-0.9747 0.9115 -0.7843 

0.1319 -0.1974 -0.1848 

2s 0.0521 -0.7674 0.0566 0.7864 0.0002 -0.7222 0.0027 

-0.6275 0.5854 0.6079 
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Table 3. Occupation numbers Xni, expansion coefRcients C"' (i=l,2,3), parameters of the s.p. wave 
functions and natural occupation numbers rini p9[ for ^^O. 

HO s.p. wave functions MHO s.p. wave functions SW s.p. wave functions 



nl a = 0.604 fm^^ a = 0.582 fm^i b = 0.731 x = 4.60 fm ijni 



0.9174 0.9576 0.9967 

Is 0.8002 -0.2686 0.9359 -0.0740 0.9940 -0.0087 0.9495 

0.2936 0.2783 -0.0808 

0.9959 -0.9928 0.9919 

Ip 0.8278 0.0008 0.8000 -0.0035 0.8370 -0.1251 0.9646 

-0.0899 0.1194 0.0214 

0.2225 0.5289 0.4743 

Id 0.1249 0.8353 0.0973 0.7562 0.0655 0.1129 0.0057 

0.5028 0.3651 0.8731 

0.1464 0.0780 -0.0315 

2s 0.0920 -0.4582 0.1774 0.9970 0.1672 0.8750 0.0057 

-0.8767 -0.0032 -0.4831 
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Table 4. Occupation numbers Xni, expansion coefRcients C"' (i=l,2,3), parameters of the s.p. wave 
functions and natural occupation numbers rjni p9[ for ^'~'Ca. 



HO s.p. wave functions SW s.p. wave functions 



nl a = 0.627 fm^^ x = 6.40 fm rj^i 



0.5323 -0.5092 

Is 0.8027 0.7972 0.9227 0.8385 0.8898 

-0.2849 -0.1939 

0.9733 0.9959 

Ip 0.8000 -0.0609 0.8718 -0.0717 0.9376 

-0.2214 -0.0560 

0.8023 0.8548 

Id 0.8005 -0.5923 0.8013 -0.5181 0.9463 

0.0746 0.0335 

-0.8440 0.8517 

2s 0.8020 0.5261 0.9558 0.4586 0.9579 

-0.1045 -0.2536 

0.2197 -0.0203 

2p 0.1005 0.5283 0.1686 -0.7750 0.0168 

0.8201 0.6317 

0.9913 0.2266 

If 0.000002 -0.1106 0.0465 0.2341 0.0127 

-0.0717 0.9454 

0.0453 0.0589 

Ig 0.1879 0.5765 0.0742 -0.2858 0.0087 

0.8158 0.9565 
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Table 5. Rms radii and mean kinetic energies per nucleon of '^He, ^^O and *°Ca calculated by using 
different s.p. wave function sets for the construction of the one-body density matrix. 

Nuclei (r2>V2 {t)/A 







(fill) 


(McV) 


A TT 


HO 


1.66 


12.33 




MHO 


1.65 


12.93 




SW 


1.67 


21.47 




EXP [26] 


1.67 






JCM [30] 


1.60 


25.86 


IbQ 


HO 


2.70 


20.08 




MHO 


2.71 


16.29 




SW 


2.60 


26.52 




EXP [33] 


2.71 






JCM [30] 


2.63 


21.34 




HO 


3.46 


38.29 




SW 


3.42 


26.94 




EXP [33] 


3.48 






JCM [30] 


3.50 


23.36 
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